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Beneath  most  overcasts,  clouds’  motions  and  rapidly  changing  optical  depths  complicate  mapping  their 
angular  distributions  of  luminance  Lv  and  visible-wavelength  radiance  L.  Fisheye  images  of  overcast 
skies  taken  with  a  radiometer-calibrated  digital  camera  provide  a  useful  new  approach  to  solving  this 
problem.  Maps  calculated  from  time-averaged  images  of  individual  overcasts  not  only  show  their  bright¬ 
ness  distributions  in  unprecedented  detail,  but  they  also  help  solve  a  long-standing  puzzle  about  where 
brightness  maxima  of  overcasts  are  actually  located.  When  combined  with  simulated  radiance  distribu¬ 
tions  from  MODTRAN4,  our  measured  radiances  also  let  us  estimate  the  gradients  of  cloud  thickness 
observed  in  some  overcasts.  ©  2008  Optical  Society  of  America 
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1.  Introduction 

Popular  and  professional  scientific  works  alike  have 
long  described  the  nonuniform  brightness  distribu¬ 
tion  of  overcast  skies  [1,2],  one  version  of  which  is 
the  CIE’s  1996  model  of  overcast  luminance  as  a  sim¬ 
ple  function  of  view-elevation  angle  [3].  Central  to 
many  such  models  is  the  assumption  that,  if  short¬ 
term  brightness  fluctuations  [4,5]  are  excluded, 
either  luminance  Lv  or  visible-wavelength  radiance 
L  is  largest  at  and  symmetric  about  the  zenith  in  op¬ 
tically  thick  overcasts  [6—9] .  Indeed,  several  mea¬ 
surement  campaigns  have  produced  data  that 
apparently  support  this  assumption  [9-13]. 

Yet  doubts  persist  about  it,  ranging  from  sugges¬ 
tions  that  thick  overcasts  have  nearly  uniform  lumi¬ 
nances  [14],  to  observations  that  some  overcasts’ 
brightness  maxima  are  in  fact  not  at  the  zenith 
[15-17],  to  an  empirical  model  that  predicts  some 
slight  zenith  asymmetry  even  for  thick  overcasts 
[18,19].  More  recently,  the  CIE  has  incorporated 
some  of  the  asymmetry  from  [18]  to  allow  for  bright¬ 
ness  maxima  around  the  Sun  in  optically  thin  over¬ 
casts  [20],  as  do  several  other  empirical  models  [21— 
23] .  Although  most  existing  data  lack  the  angular  re¬ 
solution  and  temporal  smoothing  needed  to  evaluate 


0003-6935/08/34H116-12$15.00/0 
©  2008  Optical  Society  of  America 


all  these  competing  claims  [24,25],  such  problems  are 
easily  overcome  with  digital  image  analysis  of  over¬ 
casts  photographed  with  a  fisheye  lens. 

2.  Brightness  Calibration  of  a  Skylight  Camera 

Consistent  with  our  definitions  of  “overcast”  else¬ 
where  [26,27],  here  we  analyze  brightness  distribu¬ 
tions  only  in  those  cloudy  skies  for  which  (1)  no 
clear  sky  is  visible  anywhere  and  (2)  cloud  cover  is  suf¬ 
ficiently  optically  thick  that  any  cast  shadows  are  in¬ 
distinct.  Although  this  definition  does  allow  for 
overcasts  with  thin  areas  or  a  visible  solar  disk,  nearly 
all  skies  that  we  analyzed  qualify  as  densely  overcast 
[16]  ones  in  which  the  Sun  was  indiscernible.  Our 
analysis  is  based  on  a  wide  variety  of  overcasts 
photographed  from  August  2006  to  February  2008 
at  the  U.  S.  Naval  Academy  (USNA)  in  Annapolis, 
Maryland.  From  20  different  overcasts,  we  acquired 
more  than  1300  all-sky  images  by  using  a  Nikon 
E5000  digital  camera  equipped  with  a  Nikon  FC-E8 
fisheye  converter  lens.  Although  most  of  these  over¬ 
cast  images  were  of  nonprecipitating  stratocumulus 
(Sc),  stratus  (St),  and  altostratus  (As),  we  could  also 
photograph  during  drizzle,  light  rain,  and  snow  if 
we  diligently  sheltered  and  cleaned  the  lens. 

Given  that  (1)  cloud  positions  and  optical  depths  r 
change  fairly  rapidly  during  photography  and  (2)  we 
are  chiefly  interested  in  the  angular,  rather  than  the 
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temporal,  details  of  overcast  brightness,  on  each  day 
we  took  photographs  every  30  s  for  ~30-40  min  and 
then  averaged  these  image  sequences.  At  each  pixel, 
this  averaging  reduces  brightness  fluctuations 
caused  by  short-term  changes  in  cloud  r,  fluctuations 
that  can  complicate  the  image  analysis  [28].  During 
each  photography  session,  cloud  advection  usually  re¬ 
sults  in  several  complete  exchanges  of  clouds  within 
the  camera’s  180°  field  of  view  (FOV).  Although  this 
smoothing  eliminates  the  overcast’s  cloud-level  de¬ 
tails,  it  does  preserve  persistent  features  that  are 
especially  germane  in  evaluating  the  conflicting 
claims  about  overcast  brightness  distributions. 

As  described  elsewhere,  calibrating  a  digital  cam¬ 
era  to  accurately  record  skylight  L  or  Lv  is  now  rather 
straightforward  [27,29],  as  is  determining  how  lens 
transmissivity  and  pixel  position  vary  with  zenith 
angle  for  a  given  fisheye  lens.  Note  that  because  most 
overcasts  have  a  limited  range  of  colors,  L  and  Lv  dif¬ 
fer  by  a  near-constant  scaling  factor  Lv  \L  (typically, 
the  standard  deviation  of  LV:L  <  0.5%  of  its  mean  va¬ 
lue).  Thus  the  shapes  (as  opposed  to  the  magnitudes) 
of  meridional  profiles  of  Lv  and  visible-wavelength 
L  from  400  to  700  nm  will  not  differ  appreciably  for 
a  given  overcast.  In  a  variation  on  the  spectral  calibra¬ 
tion  procedure  described  in  [27] ,  we  define  a  1  x  3  L  or 
Lv  transform  matrix  G  as 

G  =  B  pT(ppT)~\  (1) 

where  B  is  a  1  x  n  matrix  of  n  different  skylight  ra¬ 
diances  or  luminances  measured  with  a  Photo  Re¬ 
search  PR-650  radiometer  [30]  (the  training  set  of 
data),  p  is  the  corresponding  3  x  n  matrix  of  normal¬ 
ized  12bit/channel  RGB  pixel  values  for  the  same  n 
features,  and  pT  is  p’s  transpose.  Each  of  the  n  RGB 
triplets  in  p  (i.e.,  each  column)  is  multiplied  by  a  dif¬ 
ferent  scaling  factor  F  that  is  proportional  to  the 
corresponding  image’s  scene  irradiance,  and  F  is  cal¬ 
culated  from  each  image’s  photographic  exposure 
value  [31].  Our  radiance  training  set  includes  L  data 
from  five  different  overcasts,  and  it  has  n  =  232 
features  of  widely  varying  color  and  brightness.  To 
use  the  resulting  least-squares  transform  matrix 
G,  we  photograph  a  new  overcast  feature  with  the 
camera  still  operating  in  its  12bit/channel  RAW 
mode.  This  new  feature’s  reconstructed  L  is  the  only 
element  of  the  lxl  matrix  Br  given  by 

Br  =  Gpr,  (2) 

where  pT  is  the  feature’s  3x1  matrix  of  RGB  pixel 
values.  Similar  to  our  scaling  of  p’s  columns,  each 
pr  in  a  time-averaged  image  is  scaled  by  an  irradiance 
factor  F  based  on  that  image’s  mean  exposure  value. 
As  a  result  of  this  scaling,  Eq.  (2)  yields  absolute 
rather  than  relative  L  and  Lv.  According  to  Photo  Re¬ 
search,  at  specified  radiance  levels  a  properly  cali¬ 
brated  PR-650  radiometer  measures  L  and  Lv 
accurate  to  within  ±4%.  We  make  the  conservative  as¬ 
sumption  that  such  errors  follow  a  uniform  random 


distribution  (i.e.,  this  will  overestimate  many  radio¬ 
meter  L  errors).  If  this  assumption  is  applied  to  a 
set  of  camera  and  radiometer  data  not  used  to  calcu¬ 
late  G,  then  Eq.  (2)’s  root-mean-square  errors  in  abso¬ 
lute  reconstructed  L  and  Lv  are  <5.8%. 

Beyond  these  purely  radiometric  calibrations  are 
two  geometric  ones.  First  we  measure  how  light  rays 
incident  at  angles  0,  from  the  FC-E8  fisheye  lens’  op¬ 
tical  axis  are  mapped  onto  the  camera  sensor.  Knowl¬ 
edge  of  this  angular  mapping  or  projection  is  needed 
to  convert  a  pixel’s  radial  distance  from  the  center  of 
an  all-sky  image  into  the  corresponding  sky  zenith 
angle  0.  A  normalized  version  of  this  pixel  distance 
is  rn  =  rabs/r(<9;  =  90°),  where  rabs  is  the  Cartesian 
pixel  distance  from  a  target  pixel  to  the  one  at  the 
lens’  optical  axis  and  r(0i  =  90°)  is  the  pixel  distance 
from  that  axis  to  pixels  at  0i  =  90°.  Thus  0  <  rn  <  1  for 
any  pixel  in  an  all-sky  image  that  is  at  or  above  the 
astronomical  horizon. 

To  perform  this  geometric  calibration,  we  mount 
the  camera  and  fisheye  lens  on  an  optical  stage  free 
to  rotate  about  a  vertical  axis  that  intersects  the  lens’ 
nodal  point.  Next  we  position  a  point  light  source  at 
the  same  height  as  the  lens’  optical  axis  and  several 
meters  from  the  camera.  Then  we  photograph  this 
light  source  in  darkness  as  we  rotate  the  stage  in 
1°  steps.  By  relating  the  light’s  changing  pixel 
position  to  the  lens’  rotation  angle,  we  establish 
Fig.  l’s  mapping  of  The  FC-E8  lens’  slightly 

nonlinear  0i(rn)  is  close  to  the  equisolid  angle  fisheye 
projection  [sin(<9j/2)  ex  rj,  which  is  intermediate  to 
the  equidistant  (0,  ex  rn)  and  orthographic  [sinfd,)  = 
rj  projections  [32]. 


Fig.  1.  Angular  mappings  for  various  fisheye  lens  projections, 
where  9t  is  the  angle  between  a  lens’  optical  axis  and  an  incident 
light  ray,  and  rn  is  the  normalized  radius  from  the  image’s  center  to 
the  pixel  illuminated  by  that  ray.  The  FC-E8  curve  is  the  measured 
mapping  for  our  particular  model  of  Nikon  fisheye  lens. 
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Our  second  geometric  calibration  also  has  a  radio- 
metric  component.  Equisolid  angle  fisheye  lenses  are 
designed  so  that  each  sensor  pixel  images  a  solid  an¬ 
gle  of  the  same  size.  Thus  if  only  a  pixel’s  FOV  gov¬ 
erned  its  irradiance  E,  then  pixel  E  would  be 
independent  of  0,  for  an  isotropic  radiance  field.  Be¬ 
cause  the  FC-E8’s  angular  projection  is  not  quite 
equisolid,  the  solid  angle  subtended  by  (and  thus 
E  on)  its  sensor  pixels  decreases  slightly  toward 
the  image’s  edge.  Compounding  this  smaller  illumi¬ 
nation  factor  at  larger  0L  are  increased  optical  path 
lengths  within  and  external  reflectances  from  the 
lens’  glass  elements.  As  0,  increases,  the  associated 
decreases  in  transmitted  light  also  reduce  pixel  E. 

To  correct  for  these  losses,  we  map  their  combined 
effects  in  Fig.  2  as  the  lens’  effective  transmissivity 
Te {f  [also  called  its  L(di)  rolloff  [33]].  To  calculate 
TeS,  we  first  photograph  at  close  range  a  diffusely 
lit  diffuse  reflector  (a  matte-finish  photographer’s 
gray  card  outside  on  an  overcast  day)  that  spans 
much  of  the  FC-E8’s  approximately  180°  FOV  [34]. 
Because  the  card’s  reflected  L  are  quite  spatially  uni¬ 
form  for  such  illumination,  any  changes  in  image 
gray  levels  measured  radially  from  the  optical  axis 
are  due  to  the  lens  itself.  Figure  2  shows  a  smoothed, 
normalized  version  of  this  radial  profile  of  our  isotro¬ 
pic  card’s  image  brightnesses,  which  are  themselves 
proportional  to  Teff.  The  small  local  maximum  in  TeS 
at  di  ~  17.8°  [note  that  Teff(0°)  =  0.9898]  is  partly  ex¬ 
plained  by  the  combined  effects  of  the  nonabsorbing 
Fresnel  transmissivity  and  our  fisheye  lens’  illumi¬ 
nation  factors  as  functions  of  0 );  others  have  also 
measured  Te ff  maxima  at  0,  >  0°  [33].  Our  tests  re¬ 
veal  no  appreciable  changes  in  TeS(0i)  around  the 
FC-E8’s  optical  axis.  Once  we  know  an  image  pixel’s 


Fig.  2.  Normalized  effective  transmissivity  TeS  as  a  function  of  0, 
measured  for  our  FC-E8  fisheye  lens.  Teg  accounts  for  all  radiance 
losses  at  the  image  sensor  due  to  the  camera  optics,  including  il¬ 
lumination  factors,  external  and  internal  reflections,  and  absorp¬ 
tion  within  the  lens. 


di,  we  can  correct  its  brightness  by  dividing  the  ori¬ 
ginal  L  or  Lv  by  the  corresponding  value  of  Tes{Gj). 

3.  Measuring  Angular  Details  of  Overcast  Brightness 

A.  All-Sky  Maps  of  Overcast  Radiance 

We  take  two  tacks  here  in  analyzing  overcast  bright¬ 
ness  distributions:  (1)  binning  the  corrected  L  to  pro¬ 
duce  all-sky  maps  and  (2)  plotting  profiles  of 
azimuthally  averaged  L  or  Lv  as  functions  of  sky  ze¬ 
nith  angle  0.  To  explain  the  former  technique,  we  be¬ 
gin  with  Fig.  3’s  image  of  a  Sc  overcast  taken  on  4 
April  2007.  Figure  3  is  a  gray-scale  version  of  the 
average  12bit/channel  color  image  that  we  calcu¬ 
lated  from  83  individual  photographs  of  this  over¬ 
cast.  Although  we  have  made  no  geometric  or 
radiometric  corrections  to  its  pixel  gray  levels,  Fig.  3 
clearly  conveys  the  impression  that  the  circumze- 
nithal  region  is  the  brightest  part  of  the  sky. 

Figure  4  quantifies  this  impression  with  a  map  of 
Fig.  3’s  relative  radiances  Lre i  that  have  been  cor¬ 
rected  for  Teff(0i)  and  then  binned  into  5%  intervals. 
A  legend  indicates  the  gray-scale  encoding  of  Lrel, 
which  can  be  converted  to  absolute  radiances  by  mul¬ 
tiplying  it  by  Fig.  4’s  maximum  radiance  Lmax.  Fig¬ 
ure  4  and  all  subsequent  Lrel  maps  preserve 
Fig.  3’s  essentially  equisolid  projection  of  0i(rn). 
The  zenith  is  marked  by  a  black  +,  and  two  small 
black  squares  near  the  6  o’clock  position  show  the 
limits  of  the  Sun’s  position  during  photography. 
Figure  4  lists  other  relevant  details,  including  the 
photographic  session’s  time  interval  At  and  range 
of  unrefracted  Sun  elevations  /i0- 
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Fig.  3.  Gray-scale  version  of  a  time-averaged  color  image  of  a 
stratocumulus  (Sc)  overcast  photographed  at  USNA  in  Annapolis, 
Maryland  on  4  April  2007.  Averaging  of  small-scale  details  in  the 
83  individual  photographs  that  comprise  Fig.  3  makes  this  Sc  over¬ 
cast  resemble  stratus  (St). 
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Fig.  4.  Map  of  binned  relative  radiances  Lrel  calculated  from 
Fig.  3  by  using  Eqs.  (1)  and  (2)  and  corrected  for  brightness  rolloff 
in  the  original  image.  Lre i  is  Fig.  3’s  absolute  L  in  W/(m2  sr)  nor¬ 
malized  by  its  maximum  radiance  Lmax.  A  black  +  marks  the 
zenith,  and  two  small  black  squares  near  the  6  o’clock  position 
show  the  limits  of  the  Sun’s  position  during  photography.  This 
and  subsequent  maps  of  Lrel  include  details  about  the  underlying 
mean  image  and  its  corresponding  range  of  unrefracted  Sun 
elevations  h0. 


As  expected  from  [6-13],  Fig.  4’s  largest  L  occur 
near  the  zenith.  Yet  these  95th-percentile  Lre\  are 
not  exactly  centered  on  it,  nor  is  the  rest  of  the  4  April 
2007  radiance  field  especially  symmetric  (i.e.,  Lrel  va¬ 
lues  are  not  constant  along  almucantars).  Note  in 
particular  that  most  Lrel  contours  are  wider  to  the 
left  of  the  zenith  for  d  <  70°.  Because  the  Sun’s  disk 
was  never  visible  through  the  4  April  2007  overcast, 
we  know  that  its  cloud  cover  was  optically  thick. 
Thus  we  cannot  attribute  this  overcast’s  mild  asym¬ 
metry  in  Lrel  to  the  kind  of  direct-beam  scattering  si¬ 
mulated  by  the  2003  CIE  model  [20].  We  define  an 
overcast’s  asymmetry  axis  as  the  line  that  passes 
through  the  center  of  its  95th-percentile  Lre]  and 
the  zenith,  assuming  that  these  two  points  differ. 
Contrary  to  the  CIE  model’s  assumption,  Fig.  4’s 
asymmetry  axis  is  rotated  ~45°  clockwise  from  the 
Sun-zenith  axis,  the  projection  of  the  single-scatter¬ 
ing  principal  plane  onto  an  all-sky  image. 

In  fact,  Fig.  4’s  entire  radiance  distribution  is 
skewed  with  respect  to  this  axis,  a  scattering  pattern 
not  seen  in  clear  skies.  Yet  even  moderately  thick 
(say,  vertical  thickness  h  ~  0.5  km)  overcasts  make 
photons  undergo  many  more  scatterings  before 
reaching  the  surface  than  do  cloud-free  atmospheres, 
so  their  brightness  patterns  should  differ  greatly 
from  clear-sky  ones.  Multiple  scattering  in  dense 
overcasts,  especially  when  combined  with  our  tem¬ 
poral  averaging,  will  make  droplet  properties  such 
as  single-scattering  albedo  or  asymmetry  parameter 
contribute  negligibly  to  these  patterns.  Overcasts 
with  more  symmetric  Lrel  contours  such  as  Fig.  5  also 
can  have  a  distinctly  noncircular  region  of  maximum 


Fig.  5.  Map  of  binned  Lrel  calculated  from  the  time-averaged  col¬ 
or  image  [Fig.  4(a)  in  [27]  J  of  a  Sc  overcast  photographed  at  USNA 
on  30  August  2006. 


radiance  around  the  zenith.  Recall  that  all  our  maps 
of  binned  Lrei  are  averages  of  ~30  min,  so  such  asym¬ 
metries  are  far  from  momentary. 

Still  other  kinds  of  irregular  features  can  occur, 
such  as  Fig.  6’s  region  of  95th-percentile  Lrel  that 
is  neither  circular  nor  centered  on  the  zenith.  As 
in  Fig.  4,  this  feature’s  asymmetry  axis  does  not  align 
with  the  mean  Sun-zenith  axis.  In  all  maps,  the  lat¬ 
ter  axis  extends  straight  down  from  the  zenith  +, 
while  the  asymmetry  axis  in  Fig.  5  extends  from 
the  zenith  to  the  horizon  position  a.  The  routine  mis¬ 
alignment  of  these  two  axes  suggests  that  some  prop¬ 
erty  of  a  dense  overcast  other  than  its  mean  optical 
depth  governs  the  asymmetry  pattern  of  Lrel. 


USNA  Lmlx  =6.2  W/(m2  sr) 

10-5-06  9.4°  <  h()  <  15.2° 

At  =  31.3  mins 
n  =  64  images 


Fig.  6.  Map  of  binned  Lrel  calculated  from  the  time-averaged  col¬ 
or  image  of  a  Sc  overcast  photographed  at  USNA  on  5  October 


2006. 
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Fig.  7.  Map  of  binned  LiA  calculated  from  the  time-averaged  col¬ 
or  image  of  a  Sc  overcast  photographed  at  USNA  on  7  March  2007. 
Unlike  other  days  at  our  site,  snow  was  falling  then  and  covered 
the  surrounding  terrain. 


Fig.  9.  Map  of  binned  Lrel  calculated  from  the  time-averaged  col¬ 
or  image  of  a  Sc  overcast  photographed  at  USNA  on  12  February 
2008. 


As  further  evidence  of  this,  Fig.  7’s  map  of  Lrel  on  7 
March  2007  shows  an  optically  thick  overcast  whose 
asymmetry  axis  is  even  more  misaligned  with  the 
principal  plane.  On  this  day  light  snow  was  falling 
that  covered  our  site  and  surrounding  terrain  to  a 
depth  of  ~10cm.  Although  this  high-albedo  surface 
helped  reduce  the  range  of  observed  Lrel  to  ~2.14: 
1,  it  did  nothing  to  make  the  distribution  of  overcast 
Lre i  more  regular,  as  indicated  by  its  asymmetry  axis 
at  approximately  the  2  o’clock  position.  Yet  before 
making  too  much  of  snow  cover’s  effects,  note  that 
Fig.  8’s  snow-free  overcast  has  an  asymmetry  axis 
in  roughly  the  same  compass  direction.  Thus  in 
Figs.  7  and  8,  lines  drawn  from  the  zenith  to  the  same 


horizon  feature  b  are  both  close  to  their  overcasts’  re¬ 
spective  asymmetry  axes.  Other  snow-free  overcasts 
(e.g.,  Fig.  9)  have  their  regions  of  maximum  Lrel  dis¬ 
placed  just  as  far  into  the  antisolar  sky  as  does  Fig.  8. 
So  at  least  in  our  measurements,  the  distribution  ofL 
in  dense  overcasts  is  little  affected  by  large  varia¬ 
tions  in  surface  albedo,  although  these  are  important 
in  phenomena  such  as  polar  ice  blink  or  snow 
blink  [35]. 

Even  in  overcasts  where  the  Sun’s  disk  is  occasion¬ 
ally  visible,  the  region  of  maximum  Lrel  may  not  sur¬ 
round  the  Sun’s  position  or  even  be  at  its  azimuth. 
Figure  10  illustrates  the  latter  case  for  an  overcast 
where  the  Sun’s  disk  dimly  appears  in  14  of  67  photo¬ 
graphs.  Here  the  principal  plane  is  ~25°  from  the  Lrel 


Fig.  8.  Map  of  binned  Lrei  calculated  from  the  time-averaged  col¬ 
or  image  of  a  Sc  overcast  photographed  at  USNA  on  20  February 


2007. 


Fig.  10.  Map  of  binned  Lrel  calculated  from  the  time-averaged  col¬ 
or  image  of  a  Sc  overcast  photographed  at  USNA  on  13  November 


2006. 
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Fig.  11.  Gray-scale  version  of  the  time-averaged  color  image  of  a 
Sc  overcast  photographed  at  USNA  on  29  November  2006. 


asymmetry  axis,  the  simplest  explanation  for  which 
is  that  Fig.  10’s  clouds  tend  to  be  thinner  there. 
In  Fig.  ll’s  mean  gray-scale  image  of  another  Sc 
overcast,  the  Sun’s  disk  is  just  perceptible  in  an  even 
smaller  fraction  of  the  underlying  photographs  (9  of 
65).  In  Fig.  12’s  corresponding  radiance  map,  the  re¬ 
gion  of  maximum  Lrel  straddles  the  principal  plane, 
yet  it  still  does  not  surround  the  Sun  (black  squares 
near  the  6  o’clock  position).  Thus  even  in  some  thin 
overcasts,  no  simple  relationship  exists  between  the 
positions  of  Lre]  maxima  and  the  Sun. 

Table  1  summarizes  these  Lmax  details  for  19  dif¬ 
ferent  overcasts  (it  excludes  one  whose  radiances 


Fig.  12.  Map  of  binned  Lre i  calculated  from  Fig.  ll’s  color  origi¬ 
nal.  Note  that  the  brightest  area  of  this  occasionally  thin  overcast 
lies  between  the  zenith  and  Sun,  which  is  just  visible  in  ~  14%  of 
the  individual  photographs  that  comprise  Fig.  11. 


were  measured  through  a  linear  polarizer).  In  our 
maps  of  Lrel  (Figs.  4-10  and  12),  note  that  Lmax’s  ze¬ 
nith  angle  6  and  azimuth  relative  to  the  Sun  </trel  need 
not  be  at  the  center  of  the  overcast’s  95th-percentile 
radiances.  Note  too  that  <A-ei  is  measured  counter¬ 
clockwise  from  the  Sun’s  azimuth  in  these  maps  be¬ 
cause  we  are  looking  up  at  the  sky  rather  than  down 
at  the  surface.  Rows  in  Table  1  are  sorted  in  order 
of  increasing  d(Lmax),  meaning  that  these  overcasts 
exhibit  progressively  less  zenithal  symmetry  toward 
the  bottom  of  the  table. 

In  fact,  if  we  exclude  Table  l’s  last  three  rows,  con¬ 
sisting  of  overcasts  with  bright  (and  thus  optically 
thin)  horizons,  the  mean  value  of  d(Xmax)  =  11.15°, 
the  median  0(Lmax)  =  7.12°,  and  its  standard  devia¬ 
tion  s  =  9.62°.  Using  the  same  criterion  for  r/),.(,|(Xmaxj, 
its  mean  =  197.25°  and  s  =  117.93°.  A  t  test  of  the 
<9(Lmax)  departures  from  0°  indicates  that  they  are 
significant  at  the  a  =  0.016%  level.  In  addition, 
Lmax  is  nearly  independent  of  Sun  position:  the  mag¬ 
nitude  of  the  linear  correlation  coefficient  R  between 
d(Lmax)  and  mean  h0  is  <0.12.  Thus,  by  both  visual 
and  statistical  standards,  optically  thick  overcasts  do 
not  have  Lmax  at  the  zenith,  nor  do  their  Lrei  asym¬ 
metry  axes  have  any  consistent  orientation. 

B.  Azimuthally  Averaged  Profiles  of  Overcast  Brightness 
In  the  past,  many  researchers  have  depicted  overcast 
brightness  trends  by  using  meridional  profiles  of  L  or 
Lv  rather  than  all-sky  maps  [2, 11_,  16, 17,24,25] .  De¬ 
pending  on  the  technology  then  available,  such  pro¬ 
files  may  consist  of  radiometer  measurements  at 
only  a  few  6,  and  these  may  or  may  not  be  azimuth- 
ally  averaged.  Digital  imaging’s  wealth  of  high- 
resolution  angular  detail  avoids  such  limitations, 
and  it  makes  possible  entirely  new  ways  of  analyzing 
overcast  brightness.  Yet  because  meridional  profiles 
of  azimuthally  averaged  L  and  Lv  remain  a  conveni¬ 
ent  and  conventional  analytical  tool,  we  give  several 
examples  here.  Our  earlier  caveat  still  holds:  bright¬ 
ness  in  most  overcasts  is  not  constant  along  almu- 
cantars,  so  azimuthal  averaging  of  L  or  L„  tacitly 
assumes  a  symmetry  that  seldom  exists. 

Figure  H!  compares  the  daytime  Lv  profile  of  a  ty¬ 
pical  clear  sky  with  those  for  an  As  and  a  Sc  overcast; 
the  conventional  meteorological  abbreviations  CLR 
and  OVC  indicate  the  sky  state.  In  Fig.  13  and  all 
subsequent  profiles,  overcast  brightnesses  are  aver¬ 
aged  across  all  azimuths  (i.e.,  0°  <  ^>rel  <  360°).  How¬ 
ever,  Fig.  13’s  clear-sky  average  is  confined  to  two 
90°  wide  sectors  that  are  symmetric  about  the  clear- 
sky  principal  plane  (i.e.,  45°  <  <^rel  <  135°  and  225° 
<  f/vi  ^  315°);  so  its  Lv  exclude  the  solar  and  antisolar 
regions.  From  the  zenith  to  within  a  few  degrees  of 
the  horizon,  the  clear-sky  Lv  increase  monotonically 
much  as  they  would  in  a  purely  molecular  atmo¬ 
sphere  [36].  However,  Lv  decreases  for  0  >  88°,  and 
this  local  maximum  occurs  because  scattering  within 
the  near-surface  layers  decreases  direct  sunlight’s 
radiance  enough  to  offset  the  corresponding  gains 


1  December  2008  /  Vol.  47,  No.  34  /  APPLIED  OPTICS 


H121 


Table  1 .  Sun  Elevations  for  and  Positions  of  Overcast  Radiance  Maxima 


Figures 

Date 

Cloud  h0  Interval0 
Type  (deg) 

■^max 

(W /m2  sr) 

^C^max, 

(deg) 

)  4* relC^max) 

(deg) 

13,  14 

20  Aug. 
2007 

Sc 

62.7-59.7 

40.11 

1.60 

259.86 

9  Nov. 
2007 

Sc 

15.6-10.1 

5.93 

2.54 

274.58 

15 

26  Nov. 
2007 

St 

27.6-24.7 

14.9 

2.86 

247.48 

3,  4,  21 

4  Apr. 
2007 

Sc 

47.2-52.3 

39.65 

4.02 

38.69 

5,  16 

30  Aug. 
2006 

Sc 

46.4-52.1 

24.35 

4.44 

244.06 

9,  20 

12  Feb. 
2008 

Sc 

36.8-37.3 

28.52 

5.86 

176.15 

6 

5  Oct. 
2006 

Sc 

15.2-9.4 

6.19 

6.45 

334.24 

25  Oct. 
2007 

Sc 

15.9-10.7 

2.61 

6.84 

348.21 

14  Sept. 
2007 

Sc 

53.0-50.6 

33.33 

7.40 

20.55 

7,  16 

7  March 
2007 

Sc 

39.6-43.2 

24.76 

10.85 

98.84 

19  Nov. 
2007 

St 

30.1-27.6 

13.44 

11.48 

212.84 

13 

29  Jan. 
2008 

As 

31.8-32.5 

27.66 

11.76 

44.69 

10,  17 

13  Nov. 
2006 

Sc 

32.6-33.0 

37.09 

17.97 

303.43 

8 

20  Feb. 
2007 

Sc 

40.1-39.9 

22.98 

24.41 

166.31 

23  Aug. 
2007 

Sc 

30.5-37.2 

20.00 

25.42 

34.69 

11,  12, 
17 

29  Nov. 
2006 

Sc 

25.2-27.5 

39.61 

34.55 

351.40 

5  Feb. 
2008 

As 

33.9-34.5 

37.62 

69.22 

25.50 

13  Sept. 
2006 

Sc 

48.6-51.6 

28.55 

77.20 

156.93 

12  Jan. 
2007 

Sc 

27.6-29.1 

31.31 

82.94 

332.29 

“Morning  observations  list  the  smaller  h0  value  first,  whereas 
afternoon  observations  list  it  last. 


in  skylight  that  is  multiply  scattered  toward 
us  [25,37], 

In  contrast,  both  overcasts  in  Fig.  13  become  pro¬ 
gressively  darker  as  0  increases,  as  expected  from  the 
meridional  patterns  of  Figs.  3-12.  Although  one 
clear-sky  profile  of  Lv  largely  resembles  another, 
daily  changes  in  the  slopes  and  details  of  overcast 
Lv(6)  are  pronounced.  For  example,  the  ratio  of 
maximum :  minimum  L„  is  2.49:1  on  29  January 
2008,  but  it  decreases  to  2.12:1  on  20  August 
2007,  ratios  similar  to  those  predicted  by  some  of 
the  earliest  overcast  models  [2,3]. 

Yet  this  expected  feature  of  Lv  decreasing  with  in¬ 
creasing  6  is  joined  by  an  entirely  new  and  unex¬ 
pected  one:  near  the  horizon  are  paired  local 
minima  and  maxima  in  both  overcasts’  Lv,  a  common 
overcast  detail  that  is  unexplained  by  existing  radia¬ 
tive  transfer  models.  The  image  underlying  Fig.  13’s 
20  August  2007  profile  is  Fig.  14,  and  its  luminance 


2 

luminance  (cd/m  ) 

Fig.  13.  Meridional  profiles  of  skylight  luminance  L„  as  a  func¬ 
tion  of  sky  zenith  angle  8  at  USNA  for  one  clear  sky  (labeled 
CLR)  and  two  overcasts  (labeled  OVC).  Lv{8 )  in  each  profile  are 
averaged  over  all  relative  azimuths  iprei  at  fixed  6.  Although  the 
near-horizon  brightness  maximum  in  the  clear  sky  is  fairly  well 
understood  [37],  similar  extrema  in  overcasts  are  not. 

banding  near  the  horizon  seems  to  be  the  result  of 
perspective  compression  of  parallel  striations  in 
the  time-averaged  t  [27].  In  fisheye  photographs, 
these  alternating  bright  and  dark  bands  in  the  mean 
low-level  flow  resemble  barrel  staves.  Like  barrel 
staves,  the  bands  become  visually  compressed  near 


Fig.  14.  Gray-scale  version  of  the  time-averaged  color  image  of  a 
Sc  overcast  photographed  at  USNA  on  20  August  2007.  The  alter¬ 
nating  bright  and  dark  bands  are  caused  by  persistent  minima  and 
maxima  in  overcast  optical  depth  that  parallel  the  mean  low-level 
flow. 
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the  horizon  and  so  cause  the  kind  of  angularly  nar¬ 
row  maxima  and  minima  plotted  in  Fig.  13. 

Figure  15  supports  our  claim  in  Section  2  that  in  a 
given  overcast  the  meridional  profiles  of  Lv  and 
visible-wavelength  L  will  have  nearly  the  same 
shape.  This  correspondence  between  luminance 
and  radiance  will  be  closest  in  those  overcasts  with 
limited  color  gamuts,  because  the  more  dissimilar 
the  shapes  of  an  overcast’s  visible  spectra  are,  the 
larger  will  be  both  the  differences  in  its  colors  and 
Lv :  L  ratios. 

Thus  Figs.  16  and  17  show  the  same  patterns  in 
overcast  L  as  Fig.  15  does  in  Lv:  the  maximum  azi- 
muthally  averaged  L  occurs  at  the  zenith,  and  L 
steadily  decreases  to  within  a  few  degrees  of  the  hor¬ 
izon.  As  in  Fig.  13,  most  of  the  Fig.  16  and  17  over¬ 
casts  have  closely  paired  local  minima  and  maxima 
in  L  for  0  >  85°.  However,  the  time-averaged  image  of 
the  29  November  2006  overcast  (Fig.  11)  shows  none 
of  the  striations  that  readily  explain  the  near-horizon 
L  behavior  of  the  other  overcasts.  Close  scrutiny  of 
Fig.  ll’s  horizon  reveals  only  gradual  variations  in 
brightness  as  the  cause  of  Fig.  17’s  local  minimum 
at  9  =  88.25°. 

Plausible  physical  explanations  of  these  local 
minima-maxima  pairs  include  that  (1)  as  9  increases 
near  the  horizon,  subcloud  airlight  scattering  re¬ 
verses  the  trend  of  decreasing  L  caused  by  steady  in¬ 
creases  in  overcast  slant-path  r  (i.e.,  the  local 
minimum),  which  is  followed  by  (2)  greater  sur¬ 
face-based  absorption  of  this  airlight  at  the  largest 
9  where  scattering  paths  are  closest  to  the  ground 
(i.e.,  the  local  maximum).  In  scenario  (2),  multiple 
scattering  is  effectively  reduced  because  of  short  sin¬ 
gle-scattering  distances  to  the  absorbing  surface. 

2 

luminance  (cd/m  ) 
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Fig.  15.  Meridional  profiles  of  radianceL(fi)  and  luminance  Lv(d) 
averaged  across  all  <J>le\  for  the  same  St  overcast.  In  most  daytime 
overcasts,  the  ratio  Lv(6)  '.L{6)  is  nearly  constant. 


Fig.  16.  Meridional  profiles  of  L(8)  averaged  across  all  <pTe\  for  two 
overcasts  at  USNA  on  30  August  2006  and  7  March  2007.  Compare 
these  profiles  with  the  corresponding  all-sky  maps  of  Lrei  (Figs.  5 
and  7,  respectively). 

Tentative  evidence  for  this  logic  comes  from  the 
snow-covered  day  of  7  March  2007  (Fig.  16),  for  which 
there  is  no  local  maximum  in  L{9)  just  above  the  hor¬ 
izon.  Yet  the  same  holds  true  for  the  snow-free  day  of 
26  November  2007  (Fig.  15),  so  the  issue  seems  un¬ 
resolved  for  now. 

4.  Cloud  Thickness  Gradients  and  Overcast 
Brightness  Asymmetries 

What  displaces  overcast  brightness  maxima  from  the 
zenith,  a  displacement  that  seems  independent  of 


Fig.  17.  Meridional  profiles  of L(8)  averaged  across  all  <pTe\  for  two 
overcasts  at  USNA  on  13  November  and  29  November  2006.  Com¬ 
pare  these  profiles  with  the  corresponding  all-sky  maps  of  Lrel 
(Figs.  10  and  12,  respectively). 
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Fig.  18.  Schematic  cross  section  of  an  overcast  with  vertical  geo¬ 
metric  thickness  h  that  increases  nearly  linearly  in  the  direction  x. 
In  our  calculations,  we  make  this  gradient  Ah/ Ax  constant  in  mag¬ 
nitude  and  direction.  Overcast  thickness  and  curvature  are  exag¬ 
gerated  here;  6  is  the  sky  zenith  angle  for  a  surface-based  observer. 

Sun  position  (see  Table  1)?  This  is  the  fundamental 
question  implicit  in  Figs.  3-12,  and  Fig.  18  provides 
one  answer.  This  schematic  cross  section  of  an 
overcast  adds  a  new  feature  to  existing  models’  expli¬ 
cit  or  implicit  picture  of  a  uniform  cloud  deck:  a  con¬ 
stant,  unidirectional  gradient  of  cloud  thickness  h. 

The  instantaneous  distribution  of  overcast  h  and 
its  corresponding  effects  on  ground-level  irradiances 
and  radiances  can  be  complicated  indeed  [38-42]. 
However,  the  overcast  images  that  underlie  our  Lrel 
maps  (Figs.  4-10  and  12)  are  30-40  min  averages  of 
scores  of  individual  photographs,  and  so  these  maps 
are  dominated  by  only  the  most  persistent  cloud  in- 
homogeneities.  In  particular,  advection  of  the  cloud 
field  above  our  site  greatly  reduces  each  Lrel  map’s 
range  of  averaged  h  and  thus  its  horizontal  gradients 
of  h.  This  pronounced  smoothing  lets  us  realistically 
adapt  an  explicitly  horizontally  homogeneous  model 
such  as  MODTRAN 4  [43]  to  estimate  how  changes  in 
li  might  affect  the  angular  distribution  of  overcast  L. 

To  make  this  adaptation,  we  compute  tables  of 
MODTRAN’s  overcast  L  as  functions  of  h  and  0  for 
the  combination  of  solar  h0,  cloud  type,  and  cloud  base 
height  Zbase  observed  during  each  of  several  different 
overcasts  [44].  Because  MODTRAN  absolute  ra¬ 


diances  vary  negligibly  with  </>rel  for  h  >  0.5  km,  we 
can  interpolate  within  these  tables  to  estimate  how 
MODTRAN’s  predicted  Lre]  depend  on  6  and  h  at 
all  f/>re|  (recall  that  Lrel  is  absolute  L  in  W/ (m2  sr) 
scaled  by  an  image’s  maximum  radiance  Lmax).  We  as¬ 
sume  that  within  the  camera’s  FOV  (1)  zbase  is  con¬ 
stant  and  (2)  a  linear  gradient  Ah/Ax  can  be 
substituted  for  the  true  curvilinear  gradient  above 
a  spherical  earth.  However,  neither  of  these  reason¬ 
able  assumptions  is  crucial  to  our  results  below. 

Figure  19  shows  the  effects  on  Lre\  of  changing  from 
uniform  h  [Fig.  19(a)]to  a  small,  constant  Ah /Ax 
[Fig.  19(b)] .  These  maps  of  simulated  MODTRAN  ra¬ 
diances  are  based  on  Fig.  9’s  measured  h0  and  Zbase-  If 
h  is  not  spatially  uniform,  then  Fig.  19(a)’s  pattern  of 
azimuthally  symmetric  Lre\  shifts  away  from  the  ze¬ 
nith  and  develops  a  slight  asymmetry  axis  at  the 
(pTe i  where  Ah/ Ax  intersects  the  zenith  [Fig.  19(b)]. 
As  expected,  this  gradient  is  quite  small,  having  a 
value  of  ~0. 00155  or  <0.1°.  These  values  assume 
that  the  total  horizontal  distance  across  our  180° 
FOV  is  twice  the  tangent-line  distance  from  our  site 
to  Fig.  9’s  observed  zbase  =  1.4  km;  actual  horizontal 
visibility  distances  to  the  cloud  base  will  be  smaller. 
Although  we  cannot  reproduce  all  of  Fig.  9’s  details 
with  Fig.  19(b)’s  simple  unidirectional  h  gradient, 
adding  it  to  Fig.  19(a)  clearly  makes  for  a  much  more 
realistic  map. 

A  natural  corollary  to  Fig.  19(b)’s  simulated  distri¬ 
bution  of  Lre]{0, <pre\)  given  h  is  a  map  of  estimated  h 
given  the  measured  L{0,  (f>re\).  To  calculate  such  a 
map,  we  start  with  the  values  of  L(0.  </>rel)  at  each  pix¬ 
el  in  the  time-averaged  image  used  to  generate  Fig.  9. 
Then  we  interpolate  at  fixed  6  within  the  appropriate 
MODTRAN  tables  of  L(0,h)  to  find  the  value  of  h 
that  corresponds  to  each  of  these  observed  radiances. 
For  h  >  0.5  km  and  fixed  6,  MODTRAN  L  decrease 
monotonically  as  li  increases.  Thus  once  a  pixel’s  0 


Fig.  19.  MODTRAN4  simulations  of  Lrel  in  a  Sc  overcast  with  (a)  constant  h  =  0.894  km  and  (b)  constant,  unidirectional  Ah/ Ax  =  0.00155 
(~0.089°)  at  an  angle  of  34°  from  the  principal  plane.  TheLrei  scale  here  is  the  same  as  in  Fig.  9,  and  h  increases  from  0.687  to  1.101km  over 
a  horizontal  distance  that  is  twice  the  tangent-line  distance  from  our  site  to  Fig.  9’s  observed  zbase  =  1.4  km. 
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is  calculated  in  a  dense  overcast,  we  can  use  the  mod¬ 
el  L(0.  h )  to  identify  its  vertical  thickness  h  uniquely. 
Of  course,  for  an  overcast  that  actually  has  multiple 
discrete  layers,  this  h  value  will  only  be  an  equiva¬ 
lent  thickness. 

In  Fig.  20  we  estimate  the  actual  h  distribution  ex¬ 
isting  during  Fig.  9’s  measurements  on  12  February 
2008.  Not  surprisingly,  Fig.  20’s  h  pattern  is  much 
more  complicated  than  that  in  Fig.  19(b),  which  as¬ 
sumed  an  h  gradient  having  only  one  direction  and 
magnitude.  That  said,  the  measured  gradient’s  trend 
of  decreasing  h  at  (pre\  ~  130°  is  clear  in  Fig.  20  (look 
from  the  zenith  toward  location  c),  and  its  orientation 
is  reasonably  close  to  the  (p„,\  =  146°  direction  of 
-Ah/ Ax  chosen  for  Fig.  19(b).  We  estimated  the  lat¬ 
ter  direction  purely  visually  from  Fig.  9,  and  so  we 
should  not  expect  perfect  agreement  between  the 
two  gradient  orientations.  However,  much  of  the  dif¬ 
ference  arises  from  other  h  gradients  present  in 
Fig.  20,  and  the  combined  result  of  these  is  a  net  h 
gradient  that  nearly  parallels  the  principal  plane. 

Figure  21  shows  the  distribution  of  h  observed  for 
Fig.  4’s  overcast  of  4  April  2007,  and  here  the  net  h 
gradient  is  even  easier  to  see  along  an  axis  ~45° 
clockwise  from  the  principal  plane.  This  h  gradient 
manifests  itself  in  Fig.  4  as  its  Lrel  asymmetry  axis. 
Although  the  spacing  of  h  contours  in  Fig.  21  may 
seem  smaller  than  in  Fig.  20,  in  fact  the  magnitudes 
of  these  gradients  differ  relatively  little.  First,  the 
figure  legends  in  Figs.  20  and  21  show  that  their 
ranges  of  observed  h  are  nearly  equal.  Second,  the 
apparently  larger  h  gradient  on  Fig.  21’s  left  side 
is  actually  due  to  perspective  compression  of  radial 
distances  at  the  cloud  base  as  9  increases.  Clearly 
these  figures  provide  only  estimates  of  h,  and  their 
reliability  depends  on  the  accuracy  of  both  our  mea- 


Fig.  20.  Map  of  binned  overcast  thickness  h  calculated  from 
Fig.  9’s  measured  L(0,Aei)  and  MODTRAN  simulated  L(0,  h),  with 
model  parameters  set  to  match  the  conditions  observed  at  USNA 
on  12  February  2008.  The  mean  h  =  0.858  km  and  its  standard  de¬ 
viation  s(h)  =  0.0598  km. 


Fig.  21.  Map  of  binned  overcast  h  calculated  from  Fig.  4’s  mea¬ 
sured  L(9,(/>lel)  and  MODTRAN  simulated  L{9,li),  with  model 
parameters  set  to  match  the  conditions  observed  at  USNA  on  4 
April  2007.  The  mean  h  —  0.80  km  and  s(h)  =  0.0742  km. 


sured  L{6.  </;rel)  and  MODTRAN’s  simulated  L{0.  h). 
Future  research  should  incorporate  independent, 
high-resolution  measurements  of  zbase  and  h  in  order 
to  vet  our  preliminary  results.  Nonetheless,  Figs.  20 
and  21  demonstrate  the  very  real  potential  for  mak¬ 
ing  realistic  ground-based  measurements  of  overcast 
h  in  great  angular  detail. 


5.  Conclusions 

Our  research  literally  fills  in  some  angular  details 
missing  from  earlier  work  on  overcast  brightness  dis¬ 
tributions.  Short-term  temporal  smoothing  of  fisheye 
images  of  overcast  L  and  L„  lets  us  see  both  the 
small-scale  irregularities  peculiar  to  individual  over¬ 
casts  and  the  large-scale  patterns  common  to  many 
overcasts.  Naturally,  these  insights  are  possible  only 
after  Section  2’s  rigorous  radiometric  and  geometric 
calibration  of  our  all-sky  imaging  system. 

We  have  found  common  new  overcast  features  such 
as  the  asymmetries  in  L(9)  that  both  (1)  explain  some 
apparent  inconsistencies  among  earlier  papers  and 
(2)  require  new  explanations  of  their  own.  Earlier  dif¬ 
fering  opinions  about  the  angular  distribution  of 
overcast  L  and  Lv  can  now  be  seen  as  a  problem  of 
undersampling:  if  only  a  relatively  few  radiance 
measurements  are  available,  then  reconciling  their 
patterns  with  the  appealing  (but  usually  incorrect) 
assumption  of  zenithal  symmetry  is  nearly  impossi¬ 
ble.  And  although  zenith-to-horizon  profiles  of  over¬ 
cast  brightness  remain  useful,  we  now  know  a 
fundamental  caveat  about  them:  even  for  optically 
thick  overcasts,  brightness  is  seldom  constant  along 
almucantars.  So  Subsection  3.B’s  near-horizon  beha¬ 
vior  of  L(9)  notwithstanding,  meridional  profiles  of 
overcast  brightnesses  that  are  averaged  over  all 
tprA  may  obscure  important  details. 
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Thus  the  simplest  explanation  for  overcasts’  L{6) 
asymmetries  is  also  one  of  the  most  obvious:  real 
overcasts  do  not  have  spatially  constant  h,  even  when 
averaged  over  tens  of  minutes.  Instead,  changing 
gradients  of  h  appear  to  variously  tilt  and  distort 
the  time-averaged  patterns  of  L{0),  which  themselves 
depend  largely  on  paired  increases  in  6  and  slant-path 
t  through  the  overcast  [27],  If  only  this  t(0)  depen¬ 
dence  governed  overcast  brightness,  then  zenithal 
symmetry  would  indeed  be  the  norm.  This  fundamen¬ 
tal  new  observation  about  overcasts  may  well  be  use¬ 
ful  not  only  to  researchers  in  atmospheric  optics,  but 
also  to  those  in  lighting  engineering,  animal  naviga¬ 
tion,  and  photosynthetic  productivity. 
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